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Beschreibung 



von 



■•^mm^^m^ermmmmiz ^mit^em. -<S®mpu^r> aw^^a^reter 
^rw^awa***^ Daten 

-^.usuMessiaagen -von *ernmagnetischer Resonanz 

Die Erfindung betrifft einen Computer zur Auswertung von 
Daten aus Messungen von kernmagnetischer Resonanz, wobei die 
Daten wenigstens ein Relaxationssignal einer Probe enthalten. 

Die Erfindung betrifft ferner einen Kernresonanzt-dmographen 
und ein ^erfva.h>r-en ; ^uin^w von 
ker nma gne t i s cherr Res-onan z /mmhe i**wten4gstens 
Relaaea t i <^<ss|^aJ*6iS>ii^ia*^^ t^wisrd . 



ein 



Die.kelmii^e*!^ 

- NMR) wi^d'^zurein^esetzt^ein^ 

oder eine spektroskopische Information uber eine Substanz zu 
erhalten. Magnetresonanz-Bildgebung (Magnetic Resonance 
Imaging - MRI) und Magnetresonanzspektroskopie (Magnetic 
Resonance Spectroscopy - MRS) machen es moglich, eine 
regionale Hamodynamik mit Veranderungen in Blutvolumina und 
Blutzustanden sowie Anderungen des Stof f wechsels in vivo in 
Abhangigkeit von einer Gehirnaktivitat zu unte.r,s.uohen, . siehe : 
S. Pos.se«et. al.: Functdronal Magnetic -Resonance: Sfeudsbes, of 
Brain Actuation ;i:€-emi-nars.,dn Clinical- Neuropsychiatry, Vol . 
1 , No 1 , JgSta 6 . p . -7-6 -=58 8 . 

Insbesondere in der medizinischen Forschung besteht ein 
Bediirfnis, Inf ormationen uber die Gehirnaktivitat mittels 
Messungen von BlutflulS oder 

Deoxyhamoglobinkonzentrationsanderungen zu erlangen. Die 
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neuronale Aktivierung aufiert sich in einer Zunahme des 
Blutf lusses in aktivierten Gehirnarealen, wobei es zu einer 
Abnahme der Deoxyhamoglobinkonzentration kommt. 
Deoxyhamoglobin (DOH) ist ein paramagnetischer Stoff, welcher 
die Magnet f eldhomogenitat verringert und damit die 
Signalrelaxation beschleunigt. Sinkt die DOH-Konzentration 
aufgrund einer einen BlutfluB auslosenden Gehirnaktivitat, so 
wird die Signalrelaxation in den aktiven Arealen des Gehirns 
moduliert. Angeregt werden in erster Linie die Protonen des 
Wasserstoffes im Wasser. Eine Lokalisation von 
Gehirnaktivitat wird ermoglicht , indem eine Untersuchung mit 
funktionalen NMR-Methoden angewendet wird, welche das NMR- 
Signal mit einer Zeitverzogerung (Echo-Zeit) messen. Dies 
wird auch als suszeptibilitatsempf indliche Messung 
bezeichnet. Der biologische Wirkmechanismus ist in der 
Literaiiur unter dem Namen BOLD-Effekt (Blood Oxygenation 
Level Dependence - Effekt) bekannt und fuhrt bei 
suszeptibilitatsempfindlichen magnetischen Resonanzmessungen 
bei einer Feldstarke eines statischen Magnetfeldes von zum 
Beispiel 1,5 Tesla bis zu ca. 10%igen Schwankungen der 
Bildhelligkeit in aktivierten Hirnregionen. Anstelle des 
endogenen Kontrastmittels DOH konnen auch andere 
Kontrastmittel treten, die eine Anderung in der 
Suszeptibilitat hervorrufen. Durch NMR-Bildgebungsmethoden 
werden Schichten oder Volumina selektiert, die unter dem 
geeigneten Einstrahlen von Hochf requenzimpulsen und dem 
Anlegen von magnetischen Gradientenf eldern ein Mefisignal 
liefern, welches digitalisiert und in einem zwei- oder 
dreidimensionalen Feld im Meficomputer gespeichert wird. 

Aus den auf genommenen Rohdaten wird durch eine zwei- oder 
mehrdimensionale Fourier-Transformation die gewiinschte 
Bildinf ormation gewonnen (rekonstruiert ) . 
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Ein rekonstruiertes Schichtbild besteht aus Pixeln (- Picture 
element = Bildetenent ) , ,ein^o|^^| a fe^ a a±z, « s Vexeln ( = 
Volume ..etement. = Vo^tamene-iement) .....Ein Pixel 1st ein 
zweidtosionaies Bildei,emen.t, -^spiels weise ein Quadrat. 
DasBild ist aus den Pixeln zusammengesetzt . Ein Voxel ist 
ein dreidimensionales Volumenelement, beispielsweise ein 
Quader, der - mefitechnisch bedingt - keine scharfen Grenzen 
aufweist. Die Abmessungen eines Pixels liegen in der 
Groftenordnung. von W, die eines Voxels von 1mm 3 . Die 
Geometrien und Ausdehnungen konnen va-r label sein. 

Da -aus:,experimentel.len ^eaterschlrwi^j^!, von 
einer streng zweidimensionalen^Ebene ausgegangen werden kann, 
wird hauf ig-.auch hier^de^B^f f MexeZ -vexw^t, belcher 
berucksichtigt,-daBrddefcB4Jfleb:e n en^ 
dritte Dimension haben. 

Durch Vergleich des gen.ess.enen Signalverlauf s in jedem Pixel 
mit demzeitlichen Verlauf einer Modellf unktion, kann eine 
Stimulus-spezifische neuronale Aktivierung detektiert und 
raumlich lokalisiert werden. Ein Stimulus kann zum Beispiel 
ein somatosensorischer, akustischer, visueller oder 
olfaktorischer Reiz sowie eine mentale oder motorische 
Aufgabe sein. Die Modellf unlet ion, beziehungswei.se , die Modell- 
Zeit reihe , ^eschreibt die erwartete Sig-nalanderung,des 
magnetischen Re so^anz signals infolge neuronaier Aktivierung. 
Diese konnen. zum- Beispiel. mitt els^empiriseher Regeln aus 
einem Paradigma des jeweiligen Experiments abgeleitet werden. 
Wesentlich ist es, eine zeitliche Verzogerung der 
Modellfunktion gegenuber dem Paradigma zu berucksichtigen 
(trage Reaktion des Blutflusses auf neuronale Aktivierung) . 
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Es ist bereits bekannt, wie Gehirnaktivierung durch 
Aktivierungsbilder dargestellt werden kann, die aus 
kernspintomogr aphis chen Daten gewonnen wurden. Die Berechnung 
und Wiedergabe der Aktivierungsbilder ist sogar in Echtzeit 
moglich, das heiBt, ein Datensatz kann in ein Bild umgesetzt 
werden, bevor der nachste Datensatz gemessen wird. Der 
zeitliche Abstand betragt dabei etwa 1 bis 3 Sekunden. 

Eine derartige Berechnung und Wiedergabe der 
Aktivierungsbilder in Echtzeit ist in der US-Patentschrift 
5 657 758 beschrieben. Dieses Verfahren zeichnet sich dadurch 
aus, dali es eine hohe zeitliche und raumliche Auflosung 
ermoglicht . 

Ein weiteres bekanntes Verfahren ist in den Artikeln 
Jezzard, P. et al., Proc. SMPM 1993, p. 1392; Biswal, B. et 
al., MRM 34 (1995) p. 537 und Purdon, P. et al., Proc. ISMRM 
1998, p. 253, dargestellt. Bei diesem Verfahren werden ein 
Meftsignal und ein Paradigma der Messung eingesetzt. Beide 
Signale werden einer Fourier-Transformation unterzogen. 

Die bekannten Verfahren werten eine Ahnlichkeit zwischen dem 
Signal des Paradigmas und der Mefidaten aus. 

Der Erfindung liegt die Aufgabe zugrunde, ein gattungsgemafses 
Verfahren so durchzuf uhren, daB ein mbglichst hohes Kontrast- 
Rauschverhaltnis erzielt wird. 

ErfindungsgemaB wird diese Aufgabe dadurch gelost, daft ein 
gattungsgemaJier Computer so ausgestaltet wird, daB der 
Computer mit wenigstens einem Auswertemittel arbeitet, v/obei 
das Auswertemittel die Daten in wenigstens zwei Anteile 
separiert, die von einer Echozeit T E auf unterschiedliche 
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Weise abhangen. 

Die^r^n#i*ng sieht fc^^^e^Qr^^Wn^mputec zu 

X^s&Hert^wiar.d/«d^ der;*NMR- 



Signalrelaxation mit einer Zeitkonstante T 2 * = — 



an 



mehreren Zeitpunkten nach einer Anregung ermittelt. 



Bei der spektroskopischen Bildgebungsmethode handelt es sich 
vorzugsweise urn eine Echo-Planar-Bildgebungsmethode, 
insbesondere urn eine wiederholte zweidimensionale Echo- 
Biltfgefei^temel^G^e,^ Anwendung 
einer, zweidimensional-en^ best eh t . 

Sine *a*aittM^e^^ 

wiede^hdf-t^wi^rd* und^oT^gsVelatee 2 Ofc.¥i s *JT0 Otems?!?b^ragt . purch 
die«meh»i&a^he*Wd^^ 

eines Signalabf alls wird ein Verlauf des Signalabf alls in der 
Abfolge von rekonstruierten Einzelbildern dargestellt. 



Eine zweckmafiige konventionelle Echo-Planar-Methode wird als 
EPI (Echo- Planar-Imaging) bezeichnet. Eine vorteilhafte 
Implementierung der erf indungsgemaJien Methode erfolgt mit 
TURBO-PEPSI, wobei PEPSI fur Proton-Echo-Planar- 
S^eGtroseppiG-i&magi^ng- skte-ht . 



Die"Arizahl^dref^BM*de 

kodiert werden, ist abhangig von der Relaxationszeit und der 
Kodierungszeit At fur ein Einzelbild. 



Vorzugsweise wird der Computer zur Auswertung von Daten aus 
der kernmagnetischen Resonanztomographie eingesetzt, wobei 
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die Daten wenigstens ein Relaxationssignal einer Probe 
enthalten und wobei die Daten in Anteile separiert werden, 
die von einer Echozeit T E abhangen und in wenigstens einen 
anderen Bestandteil, der nicht von der Echozeit T E abhangt 
und daB die Signale, die von der Echozeit T E abhangen, als 
Aktivierungssignale erfaBt werden. 

Sine Ermittlung eines Rauschsignals ist dadurch moglich, daB 
der Computer mit wenigstens einem Auswertemittel arbeitet, 
das die Daten in wenigstens einen Anteil separiert, der von 
einer Echozeit T E abhangt und. in wenigstens einen anderen 
Bestandteil, der nicht von der Echozeit T E abhangt und. wobei 
das Auswertemittel die Signale, die von der Echozeit T E 
abhangen als Aktivierungssignale erfaBt. 

Eine Separation verschiedener Bestandteile einer zu 
untersuchenden Funktion kann durch eine Ermittlung von 
Signalen, die eine unterschiedliche Abhangigkeit von der 
Echozeit T E haben, festgestellt werden. So ist es 
beispielsweise moglich, eine Amplitude S 0 von einer 
Zeitkonstanten T 2 * und /oder von einem Rauschsignal g zu 
separieren . 

AuBerdem betrifft die Erfindung einen 
Kernresonanztomographen, der wenigstens einen 
erf indungsgemafien Computer enthalt. 

Die Erfindung sieht ferner vor, ein Verfahren zur Auswertung 
von Daten aus der kernmagnetischen Resonanztomographie, wobei 
wenigstens ein Relaxationssignal einer Probe ermittelt wire, . 
so durchzufuhren, daB die Daten in wenigstens zwei Anteile 
mit einer voneinander verschiedenen Abhangigkeit von einer 
Echozeit T E separiert werden. 
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Vorzugsweise wird das Verfahren so durchgef iihrt, daft 

JPBa'feensf^sa xgopLejags-h e1$ Echiojz e i t e n 




in*.. — wtwi- tow- -^yisiTOo^^sLvcrei^ie nit, #u^,i£a"nm'en!3Jete s 
Relaxat^^^a* eine 
^«h&irg*feg^^ der Echozeit T E erfaftt 

wird und daft das Relaxationssignal in Anteile separiert wird, 
die voneinander verschiedene Abhangigkeiten von der Echozeit 
T E haben. 

Es ist bevorzugt, das Verfahren so durchzuf uhren, daii das 
Relaxationssignal aufgeteilt wird in einen Anteil, der von 
einer*>i^^ der 
nicht von der EeMMlt T&iafcmmgt uhd, 4h Jiade r*a%teil,^er von 
der ^Echozeit -T E *afc1a.a«gt^a&^ 
wird - 

Hierbei ist ^s*be*s«lers ^m^:Cim^t,^a£^. &J m^^ns ein 
Signal ermittelt wird, das proportional zu T E exp(-T E /T 2 *) 
ist, wobei der Wert von T 2 * insbesondere durch eine 
vorzugsweise separate Fitprozedur an denselben Daten 
ermittelt wird. 

Dabei ist es besonders zweckmaftig, daft T 2 * mittels der Formel 
S = So!*exp.(-Xi*/*T2*-)^g»be^e'©htnet- wiard. 

Ferner ~i st^es >^o^t eiitelf t ^dasmvie^f^h^^ 

daft statist ische Schwankungen von AT 2 * ermittelt werden. 

Hierbei ist es besonders zweckmaftig, daft eine 
Standardabweichung o(AT 2 *) berechnet wird. 
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Es ist weiterhin vorteilhaft, daft ein Quotient <y(AT 2 *)/T 2 * 
gebildet und als Mafistab fur eine Aktivitat erfafit wird. 

Dabei ist es besonders zweckmaftig, daft eine statistische 
Abweichung einer Startintensitat S 0 ermittelt wird. 

Hierbei ist es vorteilhaft, daft eine Standardabweichung <y(£ 
berechnet wird. 

Dabei ist es bevorzugt, dafi ein Quotient ct(S 0 )/S 0 berechnet 
wird. 

Es ist besonders bevorzugt, das Verfahren so durchzuf uhren, 
daft eine statistische Schwankung von einem Rauschsignal g 
ermittelt wird. 

Dabei ist es besonders vorteilhaft, daft eine 
Standardabweichung a(g) von g gebildet wird. 

Ferner wird das Verfahren bevorzugt so durchgef uhrt , daft die 
aufgenommenen Daten in einem wenigstens zwei-dimensionalen 
Feld erfaftt werden, wobei eine Feldachse (DTE) Echozeiten T E 
erfaftt und wobei eine andere Feldachse (DTR) Wiederholungen 
von Anregungen in einem zeitlichen Abstand von T R wiedergibt 

Hierbei ist es besonders zweckmaftig, daft a(AT 2 * ) und a(g) 
durch nachfolgende Schritte bestimmt werden: 

(i) Anpassen von ilber die andere Feldachse (DTR) 

gemittelten Signalen an einen exponentiellen Zerfall 
in Abhangigkeit von der einen Feldachse (DTE) und 
Bestimmung von S 0 und T 2 *; 
(ii) Berechnung von <y(AS 0 ), a(AT 2 * ) und o(g) fur mehrere 
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Voxel und verschiedene T E mit anschliefiender Mittelung 
dieser Werte uber wenigstens eine zu untersuchende 
Region (ROI); 
( ii i )®*®tn^s®ermm3n 



o{AS) 



f \ 



■ * 
V T 2 J 



* ^ —2 "T ; 



T 2 



>0 



T* 



S 0 T 2 



r 2T E l7* 2 (o( g y 



>0 



und Ermittlung von c(AS)/S 0 als Funkti 



on von T E . 



10 Hierbei ist es besonders vo^eitoaf t,j^am^i#^# &s&n .„ von 



15 



Wei te*re «*t&e r^&sm^mK^m^n^nnd^m^t^m^e 
Weifcei*foi|!a^^ 

Unteranspru^teenvund!9dfer naeffifSi^efiden Dars^M^^feevorzugt 
AusfQhr-u.i^^s^p4#l^der--Erfi^#ung anhand von 
Beispielsrechnungen, Zeichnungen und einer Tabelle. 



er 




Von den Zeichnungen zeigt: 



Fig. 1 



Multiecho-Sequenz mit mehreren Mefisequenzen , die 
jeweils einer Spinanregung (*) folgen und rait 
-E r^assiufi g von*mefe>Be re n*.BS*r<s»ze iifee n~ T E , 



25 Fig. 2 ''^e:ine,,£r^»n.z+^ einer 

Methode des separaten Verarbeitens von Daten fur 
jede der Echozeiten, 



Fig. 3 

30 



ein experimentelles Unterschiedssignal einer 
f unktionalen Relaxationszeitanderung in einem 
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ausgewahlten Bildelement in Abhangigkeit von der 
Meflzeit nach einer Signalanregung, 

Fig. 4 AS aus verschiedenen Voxeln gemittelt iiber einige 
wenige ROIs als Funktion von T E fur 2 
representative Personen, 

Fig. 5 in einem oberen Teilbild eine Detektion von 

Hirnaktivierung in vier Schritten mittels einer 
konventionellen Bildgebungsmethode und in einem 
unteren Teilbild eine Detektion von Hirnaktivierung 
mit Hilfe eines erf indungsgemaiSen Verfahrens. 

Die Tabelle zeigt eine Zusammenf assung von experimentellen 
Beispieldaten . 

In Fig. 1 ist eine Multiecho-Sequenz mit mehreren 
Meflsequenzen , die jeweils einer -Spinanregung • (*) folgen und 
mit Erf assung von mehreren Echozeiten T E dargestellt . 

Die MeBsequenzen der Multiecho-Sequenz wurden mit der Turbo- 
FSPSI -Method? ermittelt. Jede der MeBsequenzen enthalt zwolf 
Echosignale mit Echozeiten, die zwischen 12 und 213 ms 
liegen. Die Echozeiten wurden jeweils als ein 18,3 ms 
dauerndes Zeitintervall AT E erfafit. 

Die angegebenen Werte fur die Echozeiten und die 
Zeitintervalle werden jeweils an die Geschwindigkeit der 
Datenverarbeitung angepaJit. Insbesondere konnen bei einer 
weiteren Verbesserung der Scannertechnologie die Anzahl der • 
Echosignale erhoht und die Zeitintervalle AT E verkurzt 
werden. 
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In Fig. 2 ist eine Prinzipskizze dargestellt , die zeigt, wie 
aus verschiedenen Melisequenzen jeweils ein Signal zu einer 
ersten *Echozeit be*zxdfouo*gsweise^^ oder 

Bei**der *in*F*ig. ! ' 3^da®ge^tell^:errKiirve ist- ein> Meftsignal <r(S) 
in Abharigigkeit von der Echozeit erfafit. Hierbei ist ein 
Prinzip mit einer Fitprozedur zur Aufteilung des Meflsignals 
a(S) in Beit rage, die von T 2 * abhangen und in von T E 
unabhangiges Rauschen dargestellt. Das Meflsignal cr(S) setzt 
| sich zusammen aus einern Anteil, der von einer Amplitude S 0 
abhangt, aus einen Anteil, der von einer Relaxations zei t T 2 * 
abhangt und a^^eiriem^^tfs t ar^t en^Rafus chsirg.nal^g . 

Die Er f in v dung&sTi^t^ vor ^eMie : Un^s?sthei-dung 

zwiscken«rM$%v^ ei-higm^Rab^^^durvGh' eine 

Analyse^ef:^^ 

ihrer * st^at d^jfe^s^hPe n^Wte'fet .e«M:ia|i g *z u*$ge.wi*n%i r e n . 

Eine experimentelle Uberprufung des erfindungsgemafien 
Auswerteverf ahrens erfolgt beispielsweise anhand von 
Ikernspintomographischen Untersuchungen des Gehirns von 
Versuchspersonen. In unmittelbarer Nahe des Gesichts der 
Versuchspersonen wurde eine Lichtquelle, insbesondere eine 
.Matrix v,on .Lumiaeszen,zdio.den (Light^Emi^ti.ng . J9i~ode LED) , 
pos i tdroni^er t^«and,*zu^-S i%»Rattbli t*z:en :3pge^eg t . --Dd e 
An*reg^g^f^ 8*Hz .-vEin^E^inw^^ken^der 

Signalblj^zerer^ eines 
Scanners synchronisiertes Zeitintervall von mehreren 
Sekunden, beispielsweise 5 Sekunden, an das sich ein etwa 
gleichlanges Ruheintervall anschlieftt. Bei dem Scanner 
handelt es sich urn einen Vision 1,5 Tesla Ganzkorperscanner 
der Siemens Medical Systems, Erlangen, mit 
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Magnetfeldgradienten von 25 mT/m. Ein derartiger Scanner ist 
in der Lage, Gradientenf elder innerhalb von etwa 600 ps 
umzuschalten. 

5 Als spektroskopische Bildgebungsmethode wurde TURBO-PEPSI 
(Proton-Echo-Planar-Spetroscopic-Imaging) eingesetzt. 



on: 




Eine Datenanpassung erfolgt gemafi der Exponentialf unkti 

-T IT* 
S = S Q e T E ,T 2 , 

wobei ein nicht linearer Least-square-Fit eingesetzt wird. ' 

Nachfolgend wird eine Unterscheidung zwischen Aktivierung und 
15 Rauschen mittels Multi-Echo-fMRI dargestellt. 

Die Erkennung von physiologischem Rauschen (zum Beispiel 
verursacht durch Herzschlag) erfordert ein stationares 
Frequenzspektrum, hinreichende zeitliche AuflOsung und 
20 Vorkenntnisse iiber die raumlichen und zeit lichen 
^^pharakteristika des Rauschens. Erf indungsgemaB wird eine neue 
^^Methode zur Differenzierung zwischen BOLD-bezogenen 

Variationen und anderen Fluktuationen des MR-Signals (etwa 
verursacht durch thermisches Rauschen) vorgeschlagen, die 
25 vollig ohne Vorwissen uber ein Stimulations-Paradigma 
auskommt'. Die Methode basiert auf einer Single-Shot- 
Muitiecho-Sequenz, wie der in dem Artikel von Posse, S. et 
al.. PROC. ISMRM 1998, p. 299 dargestellten Turbo-PEPSI- 
Technik. Auf diese Verof f entlichung wird vollinhaltlich Bezug 
30 genommen . 

Nach einer Signalanregung wird dessen Relaxationsverhalten in 
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Squidistanten Zeitabstanden T E aufgezeichnet . Dies wird mit 
einem zeitlichen Abstand von T R Sekunden mehrfach wiederholt. 
In^.e±ne-m -^©ken^E^^^ 

Voxels* ^n^^^eoas^na^estf e^drrnit'- dens.E-oho^eiten- T E ; in der 
.ei^en^feh^g jote ) ^^m^^^o^menm^Ta^^ndi T R in 
^de^^d'e^or Rdchtung *(*DTR) . Die Relaxation sei als 
monoexponentiell vorausgesetzt , S = S 0 exp (-T E /T 2 * J +g, mit 
einem Hardware-abhangigen Rauschen g, das wir in beiden 
Domanen, DTE und DTR, als weifi ansehen konnen. Die Werte S 0 
und T 2 * sind konstant in DTE, aber variieren in DTR: S 0 etwa 
durch Hardware-lnstabilitaten oder Blutf lufief fekte und T R 
etwa auf grund der Probanden^St i-m^afcion . ^aBdra^wen »in T 2 * 
zeigen "fcderungen «*m .l^lien^ntfftfl Wx*w*y^w;**eine 
Anderungen AS 0 und AT 2 ^ asa en. sick^^^a^ngen wie 
folgt f o>Bmul^i<e«:"eli : 



AS 



* + ~ —2 — i : — -~- 



S 0 T 2 



>0 J 



1/2 



, [1] 



wo <A> and ct(A) dem Mittelwert und der Standardabweichung 
einer GroJJe A in DTR entsprechen. Eine weitere Analyse hangt 
von der aktuellen Gr6Be der verwendeten Terme in [1] ab. Es 
ist .zwe©femaJi.i g , *&aB -f fcdie«e^er»e»elfen^^dngu^em AS 0 
sowohl- in^den ^Rufee^ais^auch, inr-dei^jct^eiw^rg^.ph»aen 
vernach>as;s.iTgbar.^ist ^aa«s^^®imeri»4m^aj^fcai»enJ;: S±n.us ) . Die 
Grofien <t(AT 2 *> und a(g) werden wie folgt bestimmt: (i) 
Anpassen des tiber die DTR gemittelten Signals an den 
monoexponentiellen Zerfall in Abhangigkeit von DTE zur 
Bestimmung von S 0 und T 2 *; (ii) Berechnung von ct (AT 2 *) und 
a(g) fur jedes Voxel und jedes T E und Mittelung dieser Werte 
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uber die interessierende Region (ROI) ; (iii, Anpassen von [1] 
rait AS 0 = 0 an diese Werte als Funktion von T E . Dies ist 
moglich, weil sich lokale Gehirnaktivierung in einem Anstieg 
von T 2 * zeigt, der eine charakteristische T E -Abhangigkeit , 
proportional zu T E e- T E /T 2 * aufweist, wogegen der Beitrag des 
weiBen Rauschens nicht von T E abhangt (s. Abbildungen) . Die 
T E -Abhangigkeit des Signals auBerhalb des Gehirns wird durch 
eine Konstante genahert. Als Validierung dieser Methode wird 
der Beitrag weiBen Rauschens mit dem Rauschen auBerhalb des 
Gehirns verglichen, unter Beachtung, da/5 a(g) aulierhalb des 
Gehirns reduziert ist. Fur eine Gaufi-Verteilung betragt 
dieser Reduktionsf aktor 0.6028. 

Visuelle Stimulationsexperimente an 4 gesunden Personen 
wurden durchgefQhrt auf einem Siemens Vision-1, 5-Tesla- 
Scanner.. Mit einer Mehrschicht-Turbo-PEPSI-Sequenz wurden 12 
EPI-Bilder (Matrixgrolie : 64 x 32 Pixel, PixelgrdBe: 3x6 
mm2) eines einzelnen FID, 90° Flipwinkel zu Echozeiten von 12 
bis 228 ms akquiriert. Eine konventionelle 
Korrelationsanalyse mittels des Sof twarepakets Stimulate 
wurde durchgefuhrt unter Verwendung eines Boxcar- 
Ref erenzvektors . 

Fig. 4 zeigt AS aus verschiedenen Voxeln gemittelt uber 
einige wenige ROIs als Funktion von T E fur 2 representative 
Personen. Die Variabilitaten aller Werte iiber ROIs waren 
klein (10-20%). Die ROIs lagen im visuellen Kortex (vc), im 
Motor-Kortex (mc) , in der weiften Substanz (wm) und auBerhalb 
des Gehirns unter Umgehung von als Geisterbildern (out) 
bezeichneten Bereichen auBerhalb des Gehirns. Die 
Filterergebnisse von [1] sind in der Tabelle zusammengef aftt , 
wo den abgekurzten ROI-Bezeichnungen die Anzahl der Voxel 
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eingeklammert folgt, ist der mittlere Korrelationskoef f izient 
tiber eine ROI, cr(g) der ROI auBerhalb des Gehirns 



r^mmmm^fS^m. s^un dr.d i e 
Eeh-ler^n«al#'erl^W.#»tem sa.nd«:f in^e.rt-»,a.ls •: eitee 



5 . S.t 



Tabelle 1 



ROI 

vc(20) - 

mc(20) 

wm(21) 

out (21) 
vc (.28) 
mc (32) 
wm (32 ) 

out (38) 



0.62 ± 0.21 
-0.11 ± 0.14 
-0.009 ± 0.19 
-0.19 ±.0.11 

0.67 ± 0. 12 
-0.22 ± 0.14 



ct(AT 2 *) /T 2 (%) 

4.3 ±0.1 
0.26 ± 0.16 
-0.001 ± 51 
•not fiitertsjed 

..3r6A±. 0 A . 1 
t 0?. 6* ± Ok 8 
r 0« 4r±4.1« 2 

^rB©t;«ear^t^d 



cr(g)/s 0 (%) 

0. 75 ± 0. 05 
0.79 ± 0. 05 
0.93 ± 0.07 
—0,. 66x.±. 0. 01 
0.42- ± 0.07 
CO . 72 ±: 0\ 06 
t^0^64i, ± 0 . 0 6 
» Of.4 5 ±'"0.01 



Bei allen^i^iso.aemvar : dej^«^rt^<>r#-a ( (.AT 2 **')^T2* "±if*den 
akt ivie^en^©«e?ln^sign^ w#g%gen-es in den 

nicht-aktivierten Voxeln zu keiner signif ikanten Abweichung 
von 0 kam. Daher ist dieser Wert deterministisch mit einer 
vernachlassigbaren stochastischen Komponente. 



Folglich kann a(AT 2 *)/T 2 * ebenso gut als Indikator fur 
regionale Gehirnaktivitat benutzt werden wie 
"•Koisr-eAat %&ms)s&&f:§i^ermen*e:^er koravent ioneM-en 
Komc&ik^t4^^m!a^^e...Zm' s >^^st)^a^.z^n - letzter, zeisgt 
<r'( AT 2 *) /T 2 ^abfer^fiur^!dren;ibe;M^^ 

Gehirnaktivitat an, so dalb eine Kenntnis eines Paradigmas 
nicht erforderlich ist. Die geringe Variability dieses 
Wertes tiber die ROIs lafct vermuten, dafi die Ergebnisse fur 
individuelle Voxel ahnlich zu den hier prasentierten sind. 
Dies erlaubt die Erstellung von o (AT 2 *) /T 2 *-Maps . Das Niveau 
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des T E -unabhangigen weiBen Rauschens liegt sehr niedrig, was 
vermuten laBt, dafi es von der Hardware herrtihrt . Das S 0 - 
Rauschen ist so klein, daB eine genauere Untersuchung des 
So-Rauschens wegen vorhandenem weiBen Rauschen schwierig ist. 

Die Erfindung sieht eine Methode zur Unterscheidung zwischen 
einer Aktivierung, insbesondere einer Hirnaktivierung, und 
Rauschen vor, wobei keine Korrelationsanalyse erforderlich 
ist. Selbstverstandlich kann die Erfindung auch in 
Kombination mit einer Korrelationsanalyse, wie beispielsweise 
einer Berechnung von Korrelationskoef f izienten, Z-Scores, 
Oder einer Anwendung des t-Tests, eingesetzt werden, um auf 
diese Weise aufgefundene Ergebnisse zu uberprufen. Eine 
Korrelationsanalyse mit zwei verschiedenen Messungen, von 
denen eine mit Stimulation und die andere ohne Stimulation 
erfolgt, ist jedoch nieht erforderlich. Eine Einbeziehung 
einer Korrelationsanalyse,. bei der Korrelationskoef f izienten 
zwischen dem Zeitverlauf der Stimulation ( „Ref erenzvektor" ) 
und Signalanderungen in Pixeln des Bildes ermittelt werden, 
kann jedoch zu Vergleichszwecken herangezogen werden. 

Hierbei aufgefundene hohe Werte des Korrelationskoef f izienten 
konnen als Aktivitatsindikator betrachtet und beispielsweise 
bei einer graphischen Wiedergabe der MeBdaten in 
Schichtbildern oder Volumenbildern als zusatzliche 
Information wiedergegeben werden. 

Die Erfindung eignet sich insbesondere fur Einsatze in 
Gebieten, in denen komplizierte Aktivierungen erfolgen. 
Deshalb ist das erf indungsgemaBe Verfahren und der 
erfindungsgemaBe Computer insbesondere zur Analyse von 
hoheren kognitiven Gehirnf unktionen wie Emotionen, Gedachtnis 
und Imagination geeignet. 
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Die Erfindung ist mit einer Vielzahl von Vorteilen verbunden. 
-eine . '<?*2n^ 

quald..tat^ven^Re-^^ ^H-i^dmrch is t es 

-m6gl!irch^~eU^ Bandbreite 
(kiirzester Kodierungszeit ) fur geringstmogliche raumliche 
Verzerrung zu verwenden und durch eine Messung einer 
optimalen Zahl von Kodierungen nach Signalanregung eine 
maximale Mefiempf indlichkeit zu erzielen. 

Das Auswerteverfahren kann in Echtzeitmessungen eingesetzt 
werden , •*^^d*ie%Re&*a«at^ 
analysieren . 




h e r a-u s g e sfc-e-xa t ,™^e£SWWam0r**ite r , *. wa s$M*tet»wte*eii*kiSft e r 
ist, eine .gewichtete Summation einzusetzen, was gegenuber 
einer Kurvenanpassung mit einer groBeren Geschwindigkeit und 
ohne Verlust an MeBempf indlichkeit erfolgen kann. Eine 
Summation, beziehungsweise eine gewichtete Summation, hat den 
Vorteil, dafi sie ein besonders robustes Auswerteverfahren 
darstellt . 

Bei al-leh^^e«Rs.U'<5h»s.pevi»s'©nen zeig-te s^ch~-ei»ne^s>ta€ke 
Aktivie-r.ung.TMm pr>im>aisen^visiael-l«eni ! fKOr.tex (Vi) und. in 
benacfebar^eraifei^a^gn .^e^gg^j^^ n<) ^^ e .^ u ^ en .^ ini mit 
TURBO-PEPSI gemessenen funktionalen Signal betragen bis zu 10 
%, abhangig von der Relaxationszeit T 2 \ der Lage und von der 
jeweiligen Versuchsperson. 

Die Anregung hat ein Maximum in der Nahe von TE = T 2 *. Bei 
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einem Vergleich von EPI und TURBO- PEPS I -Bildern mit TE = 72, 
ms wurden sehr ahnliche Aktivierungsbilder ermittelt. 

Der Empfindlichkeitsgewinn ist insbesondere fur 
Echtzeitmessungen vorteilhaft, well auch bei wenigen 
Mefiwerten eine Anderung der Relaxation wirksam ermittelt 
werden kann. Zusammenf assend kann gesagt werden, daB durch 
Multiecho-Erfassung des Unterschiedssignals eine optimale 
Empfindlichkeit bei verschiedenen Magnetf eldstarken erzielt 
wird . 

AuBerdem ist die Erfindung sowohl bei Echo-Planar-Bildgebung 
(Echo-Planar-Imaging EPI), bei phasenkodierten 
Bildgebungsmethoden als auch bei spektroskopischen 
Bildgebungsmethoden einsetzbar. 

Die dargestellten Beispiele erlautern den Computer sowie das 
Auswerteverfahren anhand von NMR-Messungen am menschlichen 
Gehirn. Selbstverstandlich konnen sowohl der Computer als 
auch der Kernresonanztomograph ebenso wie das 
Auswerteverfahren auch zur Untersuchung von anderen Proben 
lebenden Oder nicht lebenden Materials eingesetzt werden. 




Patentanspriiche 
1 . 

^Re%a»»a>i*©Kssd-gnal edmer* Probe enthalten, da- 
durch gekennzeichnet, 
daB der Computer mit wenigstens einem Auswertemittel 
arbeitet, das die Daten in wenigstens zwei Anteile mit 
einer voneinander verschiedenen Abhangigkeit von einer 
Echozeit T E separiert. 

2 . Computer -narch^n^pvctich .1 d * a srd >• u f r c *rh g * e - 
• k * e n 5n atz e , i c~h - n * e . t ^dafifdas Au-swertemit tel 
die* Daten in wenigsfceoas^ei^^ von . 
eii*8r* Etdic^t^fiiaMrtia^s^uncl .in ^n^st^ns^inen 
* andelrert*!Be1sika:h^ Te 
abtoangt undt wotoei^das ^^wertemi^^lodde^^iaXe , ...,dsLe 
von der Echozeit T E abhangen, als Aktivierungssignale 
erf afit . 

3 Kernresonanz tomograph, dadurch ge - 

k e n n z e i c h n e t, daB er wenigstens einen 
Computer nach einem der Anspruche 1 oder 2 enthalt. 

4. -Ver#ahr-en'- zur^Aui«.we , 'rt«ng- von Daten a-usgtder 

kernmagaaefeisStoer #e>Steanz , wobei^wenrgs-feens ein 
"-Re 1 ax'atinanissvrgna 1* ,e ine r Pj:&be^ri&i»t t e 1 1 , wi r d , 
dadurch gekenn zeic h- 
n e t, daB die Daten in wenigstens zwei Anteile mit 
einer voneinander verschiedenen Abhangigkeit von einer 
Echozeit T E separiert werden. 
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5. Verfahren nach Anspruch 4,dadu'rch 

gekennzeichnet,daB 
Intensitatswerte der gemessenen Daten erfaftt und in 
wenigstens zwei voneinander verschiedene Abhangigkeiten 
von der Echozeit T E separiert werden. 

6. Verfahren nach Anspruch 5,dadurch 

q e k e n n z e i c h n e t, da/5 ein AusmaB 

einer statistischen Variation der Intensitaten erfaJit 
wird. 

7. Verfahren nach Anspruch 6, d a d u r c h 

9 e k e n n z e i c h n e t, daB eine 
Standardabweichung der Intensitaten ermittelt wird. 

8. Verfahren nach einem der Anspruche 4 bis 7, d a - 

d u r c h gekennzeichn et, 
daB das Relaxationssignal in wenigstens einen Anteil, 
der von der Echozeit T E abhangt und in wenigstens einen 
Anteil, der nicht von der Echozeit T E abhangt, 
aufgeteilt wird. 

9. Verfahren nach einem der Anspruche 4 bis 8, d a - 

d u r c h gekennzeichnet, 
daB wenigstens ein Signal ermittelt wird, das 
proportional zu T E exp (-T E /T 2 *) , 1st. 

0. Verfahren nach Anspruch 9,d ad u r c h 

g e k e n n z e i c h n e t, daB T 2 * mittels der 
Formel S = S 0 exp (-T E /T 2 *) +g ermittelt wird. 

1. Verfahren nach einem oder mehreren der Anspruche 4 bis 
10, dadurch gekenn- 
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zeichnet, daB statist ische Schwankungen von 
AT 2 * ermittelt werden. 

12. .. <Vj£^.admerma^Mm®^p.i^t^2 ,rdita*d^ usr^cth 
rg*e*ic*etn «nf z e ic c«*h . n* e *jb ,«*teBi«e*ne 

Standardabweichung o-(AT 2 *) ermittelt wird. 

13. Verfahren nach Anspruch 12,dadurch 

ge kennzeichnet, daB ein Quotient 
cr(AT 2 *)/T 2 * gebildet uhd als Maiistab fur eine Aktivitat 
erf aBt wird. 

14. VerfafcvisenTnach ^^m^Mej^m^^e^rm^r.-Po^mch^ 4 bis 
13, d a -d -u r - c ~h rg t e - k * e - n *n - 

,.z M,e ^i *,c *h ^nifee^t^dSaBtfe^ 

15. Ver^ahyrenmach^A^ i 14, ta ~d. a"**d u • r -c h 
g eke n n zeich.net, daB eine 
Standardabweichung cr(AS 0 ) ermittelt wird. 

16. Verfahren nach Anspruch 15, d a d u r c h 
gekennz eich.net, daB ein Quotient 
cx(ASo)/So ermittelt wird. 

17. Verf^traenii;na^hi^ 4 bis 
-1 6 d**a «rd .„ u - r c t: h c- g*et k*e - n tr n - 
zeichnet, daB eine statistische Schwankung 
von einem Rauschsignal g ermittelt wird. 



18. 



Verfahren nach Anspruch 17,dadurch 
gekennzeichnet, dafi eine 
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Standardabweichung a(g) von g gebildet wird. 



19. 



20 



Verfahren nach einem oder mehreren der Anspruche 4 bis 
18, dadurch geke .nn- 
zeichnet, dafl die aufgenommenen Daten in 
einem wenigstens zweidimensionalen Feld erfafit werden, 
wobei eine Feldachse (DTE) Echozeiten T E erfaJSt und 
wobei eine andere Feldachse (DTR) Wiederholungen von 
Anregungen in einem zeitlichen Abstand von T R 
wiedergibt. 



Verfahren nach Anspruch 19,dadur-ch 
gekennzeichnet,daJi <x(AT 2 *> und 
a(g) durch nachfolgende Schritte bestimmt werden: 
(i) Anpassen von uber DTR gemittelten Signalen an einen 

exponentiellen Zerfall in Abhangigkeit von DTE und 

Bestimmung von S 0 und T 2 *; 
(ii). Berechnung von a(AT 2 *' und a(g) fur mehrere Voxel 

und verschiedene T E mit anschliefiender Mittelung 

dieser Werte uber wenigstens eine zu untersuchende 

Region (ROI) ; 
(iii) Anpassen von 



— < 



f \2 
T E 



2 / 



-2 



S Q T 2 



r 2T E ,7 2*JM 
< S 0 J 



und Ermittlung von ct(AS/S 0 ) als Funkt 



ion von T E . 



21. Verfahren nach Anspruch 20,dadurch 

g e k e n n z e i c h n e t, daft bei Anpa 
von or(AS)/S 0 der Ausdruck <AS 0 AT 2 *> = 0 gesetzt wird. 



ssen 
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Zusammenfassung 



5 



10 




Daten*^sMdei^^^ H ©bei die 

^a^i^^ f g^.^^i^R^^ t . i . cm§a i g ^ al , einer Probe enthalten. 

ErfindungsgemaB wird der Computer so gestaltet, dafi der 
Computer mit wenigstens einem Auswertemittel arbeitet, das 
die Daten in wenigstens zwei Anteile mit- einer voneinander 
verschiedenen Abhangigkeit von einer Echozeit T E separiert. 




Fig.3 
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